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SUMMARY 

Development of new material models for describing the "high temperature" 
constitutive behavior of real materials represents an Important area of 
research In engineering disciplines. Derivation of mathematical expressions 
(constitutive equations) which describe this high temperature material behavior 
can be quite time consuming, involved and error prone; thus intelligent appli- 
cation of symbolic systems to facilitate this tedious process can be of signif- 
icant benefit. Here a computerized procedure (SDICE) capable of efficiently 
deriving potential based constitutive models, in analytical form is presented. 
This package, running under MACSYMA, has the following features: partial dif- 

ferentiation, tensor computations, automatic grouping and labeling of common 
factors, expression substitution and simplification, back substitution of 
invariant and tensorial relations and a relational data base. Also limited 
aspects of invariant theory have been incorporated Into SDICE due to the utili- 
zation of potentials as a starting point and the desire for these potentials to 
be frame invariant (objective). Finally not only calculation of flow and/or 
evolutionary laws have been accomplished but also the determination of history 
Independent nonphysical coefficients in terms of physically measurable param- 
eters, e.g., Young's modulus, has been achieved. The uniqueness of SDICE 
resides in Its ability to manipulate expressions in a general yet predefined 
order and simplify expressions so as to limit expression growth. Results are 
displayed when applicable utilizing index notation. 


INTRODUCTION 

Development of new material models for describing the constitutive behav- 
ior of real materials represents an important area of research In engineering 
disciplines. This is evidenced by research activities, In areas associated, 
for example, with high temperature composite, reinforced concrete and geotechni- 
cal materials (ref. 3). Efforts in constitutive research involve the develop- 
ment of mathematical relationships for predicting nonlinear material response, 
derivation of material stiffness matrix appropriate for finite element calcula- 
tions (ref. 2), and finally computer implementation. The entire process re- 
quires significant manual algebraic manipulations and computer programming. 
Hence, the response time for the related efforts is quite long. As a result, 
it Is rather difficult to introduce significant changes or modifications Into 
a constitutive theory. Moreover, the outcome of the research effort may be 


affected by human errors which are often difficult to detect. In this regard, 
symbolic computation can play a major role. Immediate benefits that can be 
realized are: (1) reduced manual tedium, (2) Increased reliability of the de- 

rived equations, hence the final analysis results, (3) shortened model develop- 
ment time, and (4) Investigation of alternative functional forms. Furthermore, 
application of symbolic manipulation can provide a significant Incentive to the 
development of new constitutive theories and their finite element applications. 
However, two major obstacles arise when symbolic manipulation methods are 
applied for engineering applications; these are the number of steps In the der- 
ivation process and the problem of expression growth (refs. 1 and 4 to 7). 

Presented here Is a problem oriented software package called SDICE (Sym- 
bolic Derivation of Constitutive Equations) which Is Intended to assist In the 
derivation of potential based constitutive equations (refs. 14 and 18). The 
major features of SDICE are the automation of the equation derivation steps and 
Its ability to simplify the results so as to alleviate the expression growth 
problem. 


SYSTEM SPECIFICATIONS OF SDICE 

In deriving the material constitutive equations and matrices, six types of 
mathematical manipulations are required, l.e., 

(1) Partial differentiation 

(2) Tensor computations 

(3) Factorizations of common terms 

(4) Expression simplification 

(5) Back-substitution of Invariant and tensorlal relations 

Also, limited aspects of Invariant theory (ref. 13) has been Incorporated 
Into SDICE due to the utilization of potentials as a starting point and the 
desire for these potentials to be frame Invariant (objective). 

It has been shown that In most cases, the results obtained from direct 
application of a general purpose symbolic system, such as MACSYMA (ref. 20) are 
not useful due to the number of steps Involved In the derivation process and 
the problems associated with expression growth. For this reason, resourceful 
derivation procedures must be developed so that optimal results can be obtained 
(refs, 14 to 16). The essential features of the approach taken to address the 
above problems consists of: 

(1) A structured derivation procedure to avoid redundant steps and to min- 
imize expression growth, 

(2) Implementation of special procedures (e.g., procedures for simplifica- 
tion and pattern recognition) to facilitate the derivation process, 

(3) Expression substitution and simplification during the entire deriva- 
tion process by incorporating several levels of processing, 

(4) Automatic grouping and labeling of common factors, 

(5) Taking advantage of permutation and symmetry relationships of the 
terms Involved in each derivation step, and 

(6) As a rule-based system, intended for constitutive equation research, 
SDICE will record user defined rules and store them in a relational 
data base, whereby the Information may be retrieved, redefined and 
restored as required. 
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Some of the above procedures and techniques will be discussed, through 
application examples, in the next section. 


APPLICATION TO CONSTITUTIVE EQUATIONS FOR VISCOPLASTIC MATERIAL MODELS 

Constitutive laws provide the link between stress components on and 
strain components cij at any point in a body. These laws may be simple or 
extremely complex, depending upon the material of the body and the conditions 
to which it has been subjected. Consider the well known case of a hyperelas- 
tic material. Here, the stress and strain components are related through a 
normality structure utilizing either a strain energy or complementary energy 
function, i.e., 


aw 


ij " 3e 


13 
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e. 
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ij j 

For Inelastic material behavior the Internal state variable potential 
viewpoint is adopted, I.e., 

Q = j ,<X|j,T) 

with the generalized normality structure (refs. 8 to 10) 
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Where 0 is the complementary dissipation potential function, ch the inelas- 
tic strain, and ap the Internal state variables. Equations (4) and (5) are 
known as the flow and evolutionary equations, respectively. 

Frame Invariance (objectivity) of the resulting constitutive relations is 
Insured by requiring the potential, 0, to depend only on certain invariants and 
invariant products of its respective tensorial arguments, I.e., an Integrity 
basis. Both isotropic and transversely isotropic material symmetries have been 
considered. Transversely isotropic material symmetry Is included In the poten- 
tials of equations (1) to (3) by introducing a directional tensor d ^ d j , e.g., 

0 = (Xoj j ,ap,djdj ,T) or W = (eij.djdj). The symmetric tensor d i d j is 
formed by a self product of the unit vector d^ denoting the local preferred 
direction. 


Two viscoplastic theories have been proposed for Isotropic materials 
(refs. 11 and 12). In both theories the existence of a dissipation potential 
is assumed, and the form is taken to be, 
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Q = J f ( F ) d F + jj J g(G)dG 


( 6 ) 


where the dependence of the applied stress and internal stress (cf eq. (3)) 
enters through the scalar functions FO^j) and G(aij), respectively. 

For a material with transversely isotropic symmetry the dissipation poten- 
tial is assumed to take the form of equation (6) where now the dependence of 
the applied stress, Internal stress, and preferred direction enters through 
the scalar functions F< Ei j , d^ dj > and G(aM,didj) respectively. The stress 
dependence is now given by 


where 


F = AJ 2 + BJ 5 + Cj2] - 1 
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J 2 = 2 Vij 

J 4 = d.djEj, 

J 5 = d i d j Z jk 5: k1 
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5 4 = d.d.a^ 
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Upon application of equations (4) and (5) the resulting flow and evolutionary 
laws are: 


Gi 


f ( F) 
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where 


r SJ = ae 3 


■13 
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The computations required in the derivation of the above flow and evolu- 
tionary equations are partial differentiation and tensorial manipulation of 
Invariant relations. Direct use of a symbolic system, as our tool, in the der- 
ivation process is complicated. So, special purpose procedures have to be 
designed and Implemented in order to simplify and expedite the derivation proc- 
ess. Our work thus far has resulted in the development of several strategies, 
namely: 


(1) Store the invariant relations in a relational data base. 


(2) Implement procedures to compute tensor expressions according to the 
rules defined in tensor calculus. 


(3) Map a tensor into the domain of a scalar by utilizing a property list 
to store the variable and its subscript; thereby, allowing all differentiation 
to be treated in the same way. 

(4) Generate subscripts for Intermediate tensor variables and store them 
In the same property list in a predefined order. 

(5) Differentiation results are represented by a search tree, starting 
with the potential function as the root of the tree and its descendants with- 
out dependent relations among themselves as leaves. A separate procedure 
decides whether the function and its variables are tensorial or scalar. 

(6) Finally, check the property list and if the function Involves tensors, 
subscripts are added back for the final result according to the predetermined 
order. 

(7) Simplify the result by (a) checking the data base so as to replace any 
known invariant relation in the final expressions, (b) grouping common terms by 
factorization, and (c) identifying terms that can be written as a tensor. 

For example, when processing the Invariant J 2 the transversely Iso- 
tropic case, the definition is 


J 2-iVu 


the procedure first generates 
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following the chain-rule. Calling the procedure recursively, we have 
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and 


3 3E. . 3E,. 

EE* — U- v + — U- r 

9S k] lj 1J 3E k] L 1j 8E kl l lj 

Then checking the data base for Invariant relations and applying the tensor 
routines, we obtain from above 


2E 


1 3 


6 1k 5 jl 


Thus by combining we have, 


9J 2 

8E^ = Z 1 j S 1k S jl 

Finally, by calling the procedure Identifying the rules of tensorlal cal- 
culus and back-substituting, the result becomes E k ] . 

All the required computational rules for deriving the constitutive equa- 
tion have been Integrated Into SDICE along with rules which define tensorlal 
calculus and Invariant relations. For example. In deriving the flow law for a 
transversely Isotropic model, we have 
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From equation (6) we know that 


f 


and by tensor calculus and Indlcal notation we have 


8E 

9S 


kl 
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By Introducing constants A, B, and C Into the assumed forms of F and G 
as Indicated earlier, requires that 
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and differentiating invariant relations through the procedures we have just 
mentioned results in 


3J, 
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Next back substitution is activated, 


3J 2 3E k] ~ rt2 “kl 


3F 8J 4 


3J 4 3E kl 


2C Wl 


and 


3F 8J 5 


ai ar = B(d d,E, + d,d E , ) 

3Jg SE^i p k lp 1 q qk 

Finally, by checking the data base and substituting back all tensorial as 
well as Invariant relations defined In the table into the original expression, 
SDICE produced the following form for the flow law 

ged^ : f<ff)gk 2 (bb(3ggtl . ^ - 2gd.^jj4) + cc(6d.d^ - 2gd . j ) j j 4 + 3aa gss.^)/6gm 
where 

99tl ij : d i d i i^ ss j i i * 9!! m d n d 3 

To obtain the associated evolutionary law, equation (5) is applied, i.e., 
Instruct SDICE to take the following action 
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resulting in 


2 

gad^ : hb ged^ - g(gg) gk hb rr (bb< 3ggt2 . j - 2gd 1 jjja4) 

+ cc(6d.d. - 2gd..)jja4 + 3a a,,)/3hh 

J J 

where 

" t2 ij : d 1 d i l a ji 1 + a i 1 l d i l d j 

This agrees completely with equation (7) which was manually obtained. (We 
have adopted the convention of using double lower case letters for capital let- 
ters, and variables starting with "g" mean Greek alphabet). 


IMPLEMENTATION OF PROBLEM-ORIENTED COMPUTATION PROCEDURE 

Application of symbolic systems for material model development involves 
the symbolic solutions of systems of simultaneous equations. There are proce- 
dures under MACSYMA to solve systems of equations of different type. However, 
these general-purpose procedures when applied to our application, cause expres- 
sion growth to become uncontrollable and the results are often not useful. 

Under MACSYMA, for example, there are procedures for solving systems of 
equations, such as sol ve and 1 1 nsol ve . For certain applications, by setting 
the parameter backsubst equal to false, the results are simpler and the solu- 
tion process faster because back substitution is prevented after the system of 
equations have been triangularized. This may be necessary In very large prob- 
lems where back substitution would cause the generation of extremely large 
expressions (readers are referred to MACSYMA Manual for details). 

The systems of equations we are dealing with are primarily sparse and the 
coefficients of their variables are usually not numbers but polynomials. The 
application of Gaussian elimination which is used in 1 insol ve produces results 
that are usually lengthy and inefficient. 

A new derivation procedure is implemented into SDICE to solve the problem 
of expression growth and increase the computational efficiency. The underlying 
concept behind this Improved procedure is the identification of the smallest 
full subsystem contained within the original and then subsequent remaining sys- 
tems. Gaussian elimination is employed to solve these subsystems independently 
and sequentially Instead of the complete system. 

To clarify this procedure, consider the following sparse system of 
equations : 


2(u 2 - 
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-5u 5 x5 + 4u 4 x4 + 3u 3 x 3 + 2u 2 x2 + (u = 10 


(13) 


The procedure first searches for the smallest full subsystem In the system and 
In our case, chooses (9) and (10). The subsystem Is then solved by Gaussian 
elimination without generation of intermediate variables and we obtain: 


xl = 


V(u + t)(7ku 2 - 7ktu + 21kt 2 ) + 44u 2 
1 6ku 3 + 80ktu 2 - 88kt 2 u + 48kt 3 


x2 


(24u - 1 6t ) x 1 + 7 V(u -i- t) 
2 

1 6 u^ 


The symbols xl and x2 are then stored and treated as constants In the 
remaining equations. With the remaining 3x3 system, the smallest subsystem 
Is 1 x 1. So, an Intermediate variable zzl is generated to represent that 
part of the equation containing the previously determined variables, that Is 

zzl = 1 3u 2 x2 + 13x1 


and x3 can now be solved for Immediately 
x3 


(t 2 - 2t + Dzzl + (~8t 2 + 1 6t - 8)(kt 


u) 


1/3 


39u 

Following the same steps as before, the remaining system of unknowns is a full 
2x2 system, and two intermediate variables zz2 and zz3 are generated to 
represent previously obtained variables. 


zz2 = 


71 u 3 x3 


t 3 + 3t 2 + 3t + 1 


+ + UXl 


and 


^..3..-, (u + 2 ) x 1 

zz3 — 3u x3 + 2u x2 + ^ 


Finally, by applying Gaussian elimination, we have 
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and 


x4 


7zz3 -t- u 5 (5zz2 - 45) - 70 

9 4 

20ir - 28u 4 


4 

,. c -zz2 + 4u x4 + 9 

X 0 = y 

By comparing with the result obtained from MACSYMA , the expression size 
is smaller and the computation effort required is less, due to the generation 
of intermediate variables and the partitioning of the system Into smaller 
subsystems . 


CONCLUSION 

We have discussed the use of a symbolic computation method to automati- 
cally derive constitutive equations. The derivation steps and some of the com- 
putational methods used for potential based constitutive model research have 
been presented. It is hoped that the approach discussed here may find applica- 
tions in physics as well as other engineering disciplines. 
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